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ABSTRACT 

The virtual mass of circular cylinders immersed in a uniformly 
stratified nonhomogeneous fluid medium has been investigated both theo- 
retically and experimentally. The experimental equipment consisted of 
a reservoir filled with one inch thick layers of salt water solution. 
The density uniformly decreased with elevation (stable stratification). 
The fluid motion about the cylinder was provided by the standing waves 
created with a plunger. The resulting velocities and accelerations 
about the cylinder were harmonic. The theoretical analysis consisted 
of the calculation of the inertia coefficient through the evaluation of 
the total kinetic energy of the fldidPediun. Use was made of the psuedo- 
potential flow solution for the stratified flow about a circular cylinder. 
The results have shown that the inertia coefficient for circular cylinders 
immersed in a uniformly stratified media is about 25% larger than that for 


a cylinder immersed in an homogeneous medium. 
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hag Introduction. 

The purpose of the present investigation is the determination of 
the virtual mass coefficient for circular cylinders immersed in a uni- 
formly stratified nonhomogeneous medium. Virtual mass, which is a 
combination of the mass of the body and the added mass resulting from 
the acceleration of the fluid particles, is of considerable importance 
in the determination of the resistance of bodies accelerating from rest 
or decelerating to a full stop. It is for this reason that numerous 
studies were conducted to elucidate the concept of virtual mass and to 
determine the appropriate virtual mass coefficients for various shapes 
of bodies in homogeneous fluid media. 

It is necessary to briefly describe various definitions of the 
added mass concept prior to a description of the work undertaken here- 
with. As a body is accelerated from rest in a fluid otherwise at rest, 
the fluid medium surrounding the body is also accelerated. The signif- 
icant difference between the acceleration of the body and that of the 
fluid particles is that while the solid body has constant acceleration, 
the fluid particles have accelerations varying from that of the solid 
body to zero at an infinite distance from the body. In order to deter- 
mine the acceleration of the solid body for a given force F exerted on 
the body, it is, therefore, necessary to know not only the mass of the 
body but also the integral contribution of the fluid medium surrounding 
the body. If a volume integral is taken of the local accelerations 
times the infinitesimal fluid masses and expressed in terms of a lump 
sum of fluid mass times the constant acceleration of the body, then that 


lump sum of the f luid mass may be regarded as the mass attached to the 
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solid body or as increasing the mass of the solid body by this additional 
mass of fluid. It is for this reason that the mass contributed to the 
solid body is called the added mass and the sum of the two masses is 
called the virtual or apparent mass. In other words, the virtual mass 

is the quotient of the force required to produce the accelerations 
throughout the fluid divided by the acceleration of the body. In 


mathematical terms, Newton's second law may thus be written as 


F=(Mt+ma 
in which M is the mass of the solid body, m the added mass, and &% is 
the acceleration of the solid body. Obviously, in the case when the 
body is immersed in a fluid medium of very small density such as air, m 
is negligibly small and thus the equation of motion may be written in 
its customary form of F =Mc& . In the case where the body is immersed 
if a fluid of relatively large density such as water, then m becomes 
quite aignreieane and its neglect may lead to extremely large errors 
depending on the actual mass of the body. As an example, it is easy 
to verify that the acceleration of an air bubble will come out to be 
infinitly large if the added mass is neglected, whereas its true value 
is twice the acceleration of gravity as would be found through the in- 
corporation of the added mass into the equation of motion. 

In operations in large liquid bodies such as in the oceans, one 
must take into consideration not only the contribution of the added mass 
in determining the forces acting on the body, but also the effect of the 
variations in density of the fluid medium on the added mass itself. 

The studies carried out on added mass in the past were confined to homo- 
geneous fluid media. The advent of underwater technology and the recent 
progress made in understanding the motion of stratified fluids made it 


clear that the engineer must have a better understanding of the virtual 
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mass of various shapes of bodies in nonhomogeneous fluids. It is 
apparent that there are an infinitely large number of combinations of 
shapes of bodies and types of liquid nonhomogeneities to be considered. 
For a scientific approach to the problem, however, one must consider 
the experimentally and mathematically amenable cases such as the motion 
of a right circular cylinder in a uniformly stratified fluid medium. 

It is for this reason that the project described herein investigates the 
virtual mass of three sizes of circular cylinders in two different but 
uniformly stratified salt-water solutions through the use of standing 
waves. The results show, as anticipated theoretically, that the vir- 
tual mass coefficients in stratified fluids are about 25% larger than 


those in homogeneous fluids. 
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II. Previous Theoretical and Experimental Studies. 

The concept of virtual mass was probably first recognized by DuBuat 
in 1786 while studying the sources of errors in the nOCtonmed pendulums. 
Bessel, in 1826, also recognized the effect in his studies of the same 
problem. Analytical solutions to the problem were found by Poisson and 
Green in 1832 and 1836. Froude, in his studies of ship oscillations in 
1861, also noticed the added mass effect and took it into account in his 
inertia coefficient. 

Since 1920, there have been a number of papers and texts published 
concerning added mass, virtual mass, or inertia coefficients. Lamb (1)? 
developed some of the concepts and mathematical background of the virtual 
mass effect. Landweber and Yih (7) compared the solutions of Lagally, 
Taylor, and Kirchoff showing their correlation and simplying their proofs. 
Texts by Rouse (3) and Robertson (2) give clear definitions and theoreti- 
cal developments of the phenomenon. 

Much experimental work on virtual mass has been done in the past 
15 years. Keulegan and Carpenter (5) determined inertia coefficients 
for circular cylinders and flat plates through the use of standing 
waves. McNown (9) carried out similar studies for lenticular cylinders 
and flat plates. Sarpkaya (13) conducted studies on the added mass of 
lenses and parallel plates. Brater et. ai. (8) dealt with the determina- 
tion of the virtual mass coefficients of barges of different geometries. 
All of these studies have tended to support the theoretically determined 


virtual mass coefficients for low velocity and small amplitude motions. 
lL etilees in parenthesis refer to bibliography Section VIII. 
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As the velocity and amplitudes were increased, departures from the 
theoretical values were observed. The theoretical analyses invaribly 
assumed irrotational flow in an ideal, homogeneous, incompressible 
fluid. Thus it was to be expected that experimental values would 
differ from those obtained theoretically as separation, vortex forma- 
tion, and turbulence accompanied the flow following the initial periods 
of motion. Sarpkaya (13,14) carried out extensive theoretical and ex- 
perimental studies of the virtual mass and time-dependent resistance 
for flows with vortex separation. 

The previous experimental studies carried out with homogeneous in- 
compressible fluids with bodies either undergoing oscillations of small 
amplitude or accelerating from rest have shown that the theories of 
Lagally, Taylor, and Kirchoff are substantially correct provided that 
the fluid motion does not continue long enough in one direction to 
cause separation. 

The work on the corresponding problem in non-homogeneous fluids, 
however, has not been as prevalent. Yih (4,12) and Long (10,11) ap- 
pear to have done most of the work in this field in the past ten years. 
The most important discovery with regard to virtual mass effect is 
the phenomenon of "blocking" discussed by Yih (4). 

The experiments of Long (11) and Yih (12) indicated that for 
small amplitude motion stagnation zones occur upstream and sometimes 
even downstream from a barrier or a moving body. This phenomenon of 
stratified flow would seem to manifest itself in the form of a "frozen" 
volume of fluid following the motion of a solid body. Theoretically, 
there appears to occur a critical regime in which blocking takes place. 
If the stratification is slight, or the motion rapid, the effect is 


not noticed. In the present study, the blocking effect per.se. was 
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not studied, however the relationship between this effect and the virtual 
mass effect was necessarily taken into account. 

From the foregoing, it should be apparent that the blocking would 
increase the virtual mass coefficient. The validity of this argument 


has been shown through the experiments presented herein. 
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III. Experimental Equipment and Procedurea. 
Equipment 

For the experimental study, a one inch thick plexiglas tank with 
outside dimensions 24" by 24" by 36'' was conatructed. A * ineh hole 
was drilled and tapped for pipe threads in the center rear of the tank 
bottom. A three inch diameter circular aluminum flat plate was inatall- 
ed with a 3/16" clearance over this hole. The purpose of the hole was 
to permit slow filling of the salt-water in layera of varying denaity 
from the tank bottom. 

A pressure tap was drilled in the side of the tank 4 inechea from 
the top and centered 12 inches from either aide of the tank. This tap 
was made for use with a differential pressure transducer in recording of 
wave heights. A mounting bracket for the tranaducer was installed on 
the outside of the tank in a position such that the zero level of the 
transducer was located at the level of the pressure tap. One side of 
the transducer was connected to the tank tap and to a length of tygon 
tubing in the shape of a standpipe. The standpipe was used ag a cali- 
bration device. The other side of the differential transducer was left 
open to atmosphere for the purpose of recording gage pressures, 

A heavy steel box frame was constructed around the tank and its 
supporting table. The base of the frame reated on the concrete floor of 
the building. The top of the frame was securely bolted to a concrete 
pillar, The frame was designed so that the height of the cylinder arm 
could be adjusted in two inch increments, This was accomplished by means 
of an aluminum bracket which bolted to both sides of the steel frame, 
The frame was drilled and tapped at two inch intervals from the level of 


the top of the tank to the top of the frame, (See Fig. 2), 
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Figure 2. A schematic drawing of the test apparatus 
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The cylinder arm was attached to the bracket by a gimbal bearing 
with only one degree of freedom. The upper extension of the arm was 
clamped to the steel frame with a "C" clamp. The cylinders themselves 
were attached to the arm by another aluminum frame. This frame consist- 
ed of two 1" by 1/8' streamlined aluminum arms, 20"' long. The arms were 
bolted to a top piece which consisted of 1" by %"' aluminum bars bolted 
together in the form of a channel beam. The top piece was 21.5" long 
as were the cylinders. Thus, when the cylinders were bolted to the 
arms, a 20" by 21.5" rectangular rigid form was produced. This was in 
turn bolted to the cylinder arm at the center of the top piece. 

The cylinder arm was a 3/4" by 7/8" steel rod. A one inch test sec- 
tion was notched into the arm at a point six inches above the attach- 
ment point of the cylinder frame. The cross section of the test section 
was 3/8" by 7/8". A four gage bridge was formed at this section using 
Microsensor strain gages with a gage factor of 150. The bridge was 
completed at the test section and the leads brought out in a shielded 
cable to one channel of a two channel Sanborn amplifier/recorder system. 

Calibration of the bridge was accomplished through a pulley, weight, 
and string system. One pulley was attached to a heavy brass block and 
set into the empty tank. Another pulley was clamped to the aluminum 
cross-bracket. A string was attached to the middle of the cylinder 
and passed over both pulleys. A weight pan was built to attach cali- 
bration weights to the string. Thus a horizontal force was produced on 
the cylinder equal in magnitude to the weight of the pan and weights. 
(See Fig. 3). 

The tank itself was mounted on a heavy 4" by 4° table with a hole in 


the center rear of the table. A six inch long and %-inch diameter pipe 
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was attached to the hole in the bottom of the tank, and supported by a 
pipe bracket bolted to the table. This pipe extended below the table 
and was fitted with a valve and tubing fittings for filling or draining 
the tank. (See Fig. 2). This valve is shown as valve A. The filling 
tubing ran to valve B located at the front of the table. Tubing from 

B led to a tee connection at D which was the supply spigot from the 
building water main. From spigot D tubing was led through another tee 
either to the filling bucket or to a washdown line. Valves C and E were 
isolation valves for their respective lines. 

A simple filling and mixing bucket was made by punching a hole in 
the bottom of an old five gallon kerosene can which had been thoroughly 
cleaned. A fitting was attached to the bottom of the bucket and tubing 
run from there to the filling system as shown on Fig. 2. Thus, as can 
be seen from the diagram, by properly manipulating the valve arrange- 
ments, the tank could be filled either from the building supply with 
fresh water, or from the bucket with salt-water. 

Salt-water solutions were made in the bucket by filling the bucket 
with fresh water to a height corresponding to a one inch layer in the 
tank, then adding measured quantities of salt. In order to increase the 
pressure head while filling from the bucket, the bucket was raised about 
six feet above the level of the table. The salt solution could then 
drain by gravity into the bottom of the tank. 

A standing wave was created in the tank by means of a 20 inch by 8 
inch piece of % inch plywood. Using this as a paddle of sorts, a wave 
was gently started by hand. The amplitude of the wave was increased by 
applying a slight lift or push to the paddle in phase with the natural 


motion of the standing wave in the tank. The plywood was left floating 
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in the tank to help damp out any secondary waves, and because lifting the 
paddle out of the tank tended to create secondary waves disturbing the 
basic standing wave pattern. 

All readings were taken on a two channel Sanborn amplifier/recorder. 
This instrument consisted of a Hewlett/Packard model 7712-B recorder 
with Sanborn model 350-1100C preamplifiers. 

Procedures 

Prior to conducting the experiments a careful calibration procedure 
was followed. The first item requiring calibration while the tank was 
empty was the strain gage bridge and recorder system. This calibration 
was accomplished by the pulley system previously described and shown in 
figure 3. With no weights or weight pan on the string, and all tension in 
the string relaxed, a careful zero balance was obtained on the recorder. 
Then an accurate weight of 0.1 pounds, including the weight of the pan, 
was suspended from the string. The needle deflection on the recorder 
was adjusted by means of the preamplifier gain control so that at an 
attenuation setting of 10, a deflection of 25 millimeters (half of full 
scale) was indicated. The weight was then removed and the system zero- 
balanced again. The above procedure was followed again and checked to 
insure that the deflection for the same weight and attenuation setting 
remained constant. All higher attenuation settings were then checked to 
insure linearity. When this had been accomplished for 0.1 pounds, the 
weight was changed to 0.2 pounds and the procedure followed again. Once 
the gain was set during the initial calibration, it was locked in place. 
Calibration checks were made on both sides of the cylinder, that is the 
pulley system was shifted to the other end of the tank and calibration 
repeated for force in that direction. Calibration remained the same for 


either side of the cylinder. 
25 


The tank was then filled to a depth of 19 inches with fresh water. 
This allowed for one inch of water above the pressure tap in the side 
of the tank. The transducer pressure side was filled with water from 
the tank and all air pockets removed from the tubing and transducer 
cavity. The transducer output was connected to the other channel of 
the recorder. The recorder was again zero balanced with both the tank 
valve and the standpipe valve open. The tank stop was then shut and a 
water column of one inch above tank level was impressed on the trans- 
ducer by means of the tygon tubing standpipe. The amplifier gain was 
adjusted for half of full scale deflection with the attenuation setting 
at 2. The gain on that channel was locked and checks of linearity were 
made similar to those conducted for the strain-gage arm. Water heights 
of up to one inch above and below the level of the tank were tested, 
and the calibration was found to be linear. 

Following the completion of the calibrations, a complete system 
check of the experimental apparatus was made. With fresh water in the 
tank, runs were made with the 1.0, 1.5, and 2.0 inch diameter cylinders 
to see how closely the experimentally determined homogeneous fluid 
virtual mass coefficients agreed with those predicted theoretically. 

To do this, a standing wave was generated as previously described. Once 
the fundamental standing wave with no harmonics had been set up, the force 
acting on the cylinder and the wave heights at the end of the tank were 
simultaneously recorded. 

Immediately after each experiment, the transducer calibration was 
checked. It was found that the gain setting needed no changes through- 
out the series of experiments. The water was then drained out of the 


tank and the strain-gage arm calibration checked again as above. Again, 
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no change was required from the initial calibration setting in spite of 
changing of the cylinders. 

When all the runs had been completed with fresh water, the water was 
drained out of the reservoir, calibrations were repeated, and the strati- 
fied layers were introduced into the tank. This was accomplished by 
slowly filling the water-salt solution from the bottom of the tank through 
the filling system previously described. The first layer to be intro- 
duced was two inch layer of fresh water. The solutions for the second 
and succeeding layers were prepared by filling the bucket and adding a 
measured quantity of salt (10 milliliters or 30 milliliters) and mixing 
thoroughly. The quantity of salt was increased by equal amounts for each 
succeeding layer, that is, 10 ml. for the second layer, 20 ml. for the 
third, 30 ml. for the fourth, etc. As each solution was mixed in the buck- 
et it was dyed with red, blue, or black ink, or on occasions was left 
clear for contrast. The purpose of the dye was merely to make the layers 
visible. Once the salt, water, and dye had been mixed and the salt com- 
pletely dissolved, the specific gravity of the solution was determined 
with a hydrometer. The bucket was then raised to an elevation 6 feet a- 
bove the level of the reservoir, and the valves opened to allow draining 
of the bucket into the bottom of the tank. 

The flat plate over the filling hole allowed the slow radial diffu- 
sion of the entering solution, preventing premature or undesirable mix- 
ing of the successive layers. As the new, more dense layer entered the 
bottom of the tank and spread radially from the filling hole, it lifted 
the less dense layers above it. This process was continued until the 
tank had been filled to a water depth of approximately 19 inches. Once 


the filling process was completed, the tank was allowed to sit for a day 
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Figure 4(a) Stratified Layers 





Or so prior to the start of the tests. This waiting period allowed mole- 
cular diffusion to take place producing a uniform density gradient rather 
than a step layer effect. Each one inch layer took approximately 20 
minutes to prepare, and an hour to drain into the tank. The entire fill- 
ing procedure took a period of three days. 

Once the tank was filled with the stratified salt-water, the cali- 
bration of the transducer was checked as previously described. A stand- 
ing wave was set up in the tank. The forces acting on the cylinder and 
the wave amplitudes were recorded. Two density gradients and three 
different cylinder diameters were used in the experiments. From these 
data the virtual mass coefficients of cylinders in both the homogeneous 


and stratified fluid media were calculated. 


Zo 


IV. Theoretical Considerations. 
1. Added Mass and Inertia Coefficients (Homogeneous Fluid) 
The force exerted on a circular cylinder by an inviscid fluid 
of velocity U and acceleration dU/dt may be determined as follows: 
Letting U be the undisturbed velocity of the ambient flow, 
the velocity components for the flow about the cylinder are given by 


(See Lamb) 


u =U (kaa ayo) 


it 
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The pressure at any point in the fluid medium reduces to: 
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Writing the equation of momentum in its familiar form 


pe - 9% dA + § | (up-ub dan + \ (anpaar (3) 


=~ 6 
and noting that for the case under consideration, u,* Uy at infinity, 


i.¢., the flow is symmetrical, we have 


a = ~$h% JA + (P= P2) dx (4) 


It is apparent that the force acting on the cylinder in the 
absence of a pressure gradient is given by the first term on the right 
hand side of equation (4). This corresponds to the case where the fluid 
ig at rest and the cylinder is moving with the time dependent velocity U,. 
Because of this obvious difference between the acceleration of flow about 
a cylinder at rest and the acceleration of the cylinder in a fluid at 


rest, two types of force coefficients are employed in the literature. 
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The normalized force coefficient for the case corresponding to the ac- 
celeration of the fluid is called the inertia coefficient, and the one 
corresponding to the acceleration of the cylinder, the added mass co- 
efficient. As will be seen later, the two coefficients differ only and 
always by a constant equal to 1.0. 

Inserting equations (1) and (2) into equation (4) and performing the 


indicated integrations, we have 
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Normalizing equation (5) by "displaced mass times acceleration", 


yields 
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where Cums is known as the inertia coefficient. It is obvious that 
in the case of a circular cylinder accelerating in a fluid otherwise 
at rest k* = 1, k* being the added mass coefficient. For other shapes 
of bodies, k* differs from unity and Cc reduces to: Cs =k* +1. In 
passing it should be noted that for a sphere CT = 1.50, provided that 
the flow about the sphere is unseparated. 

In the present investigation, the cylinder is at rest and the theo- 
retical value of the inertia coefficient for an homogeneous working 
fluid is Ce = 2.0. This will later be compared with those obtained 
experimentally and it will be shown that there is excellent agreement 


between the theoretical prediction and the measured values. 


2. Added Mass Coefficients (Nonhomogeneous Fluid). 


The equations of motion for a nonhomogeneous fluid will first 
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be summarized prior to the evaluation of the added mass coefficient for 


a circular cylinder immersed in a uniformly stratified fluid medium, 


The equation of continuity may be written as 


/ / 
ou 4 ow 0 (7) 
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/ 
The stream function yv is given by 
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The equations of motion, in terms of u' and w', become 
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Writing the Bernoulli constant Hy as 
2 
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equations (10) and (11) may be written as 
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The foregoing equation was derived by Long (11) and Yih (4). 

Equation (14) may be used in the study of a class of psuedo- 
potential flows. Ome such flow case is the motion of a circular cylinder 
in a uniformly stratified fluid. The stream function for the case under 


consideration has been given by Yih (4) as follows: 
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in which Ste r/Tos ro being the radius of the cylinder, C, = U, the 
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velocity of the cylinder, and 
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The added mass coefficient may be calculated through the use 
of equation (15) by evaluating the velocity components and performing 


the integration for the kinetic energy as 
~e oT 
a 
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| = 2 age (17) 
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Y\ is the added mass sought. 


The velocity components in polar coordinates are obtained from 
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equation (15) as 
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(21) 


Obviously, the integration of the terms involving sin © sin 
30, cos 9 cos 30 through the use of equation (17) between the limits 
indicated shows that these terms have no contribution. Furthermore, the 
integration of the terms involving cok® @ and ain’ @ cancel each other 
out due to the opposite signs of their coefficients. The remaining terms, 
i.e., the first, second, and the last term on the right hand side of 


equation (21) yield 
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ee PS ae 
2 a | ° 
Teste am +5 (Tho 42 athe) (2) 


Hence the added mass YY\ reduces to 


6 
24 TPO? fe 
mM = T fot + ae Ce Ue (23) 


Inserting the value of on. from equation (16) into equation (23) and 
ya 
dividing by the displaced mass of the cylinder i.e., TY Fo Yo we 


have the added mass coefficient k*, as 


 * 1+¢(1-22)* (24) 
o daz 
It should be noted that for a density stratification where 
‘ig is decreasing with elevation (stable stratification), J P/4 > 
is negative. 
The inertia coefficient oy » which was given by Cm = ke +1, 


now reduces to 
= n -& 28)" 
Cy oy fo dz (25) 


The results obtained from equation (25) will be compared in the next 
section with those obtained experimentally. It suffices to note that. 
in stratified flows the inertia coefficients depend, unlike homo- 
geneous fluids, on the density of fluid, the radius of the cylinder, 


and of course on the density gradient. 


BD 


on Characteristics of the Experimental Waves. 

The region under the nodal area of a standing wave that may be 
realized in a rectangular tank provides a velocity field of simple 
harmonic motion in the velocity component U. 

Taking the x-axis in the plane surface of the undisturbed water, 
the 2 axis vertical and upwards, and the origin at one end of the tank, 
the surface elevation, measured from the undisturbed level, is given 


by Miche (6) to a second order of approximation in the amplitude as: 


| a*k atk 
lh = 0 Cos k x $1 at = 4 N, cos ZX = “4 Nie Cos 2KX Cos z rt (26) 


cosh 2K 
where ON, = Gin 2k 


and a cosh kH (cosh Z2kH + 2) 
No ‘ enh? KH sinh k 


where k = T/ » L being the length of the basin; Os 20 /- » IT being 
the period of oscillation; H the depth of water; and a the semi-wave 
height, that is the mean value of the extreme end deflections in a cycle. 


The period is related to k and H by the expression 
g* = g k tanh KH (27) 
The surface elevation at x = 0, at any time, becomes 
. ale = q? k cos 2et 
h = asinot + “a Ny A Nz Cos (28) 


The expressions for the particle velocities, within the order of 


the approximations considered, are given by Miche (6) as 
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~92% Cosh k (2+H) 


sin kx ¢ 
ae bh ik wi ba os ot 


~3 ga*k* gosh 2k (#+H) .. | 
4S sinh® RH emh2kH nex sin Zot (28) 


ae ee 
ov Cosh KH 


+2 gatk* cinh 2k (2+H tice kX Bin cele Co) 
+ we See hie = 7 


At the vertical plane through the midsection of the tank, i.e., at the 


Qos KX tos at 


plane x = L/2 or kx = TY /2, the velocities become 





amie ak Gosh k (2+H) i 
UU —- oth let Com (31) 


and 


2 ga. S sinh 2k (2+H) 
Wa = a “2 


on > ae. ae Zet 
;” satetenbaee - 


Obviously, the horizontal component of velocity at the channel mid- 
section is simple harmonic. The vertical component is also simple har- 
monic except that the frequency is twice as large. The effect of the 
vertical velocity is negligibly small and decreases with both the wave 
height and the lowering of the cylinder in the reservoir. The majority 
of the experiments were made with’a” less than 0.5 inch. For these cases 
the ratio of the amplitudes of vertical and horizontal components of 
velocity, Wm sven , at the cylinder location was found to be 0.004. 

The acceleration at the nodal plane may be obtained from equation 
(31) 
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UYU. gak cosh k(e+H) .. 
9 Ane ee wae Sin st (33) 


cu 


Clearly, acceleration reaches its maximum value when the velocity 

reaches its minimum value of zero and vice-versa. This observation con- 
stitutes one of the most significant points of the experimental proce- | 
dure. It is a well known fact that the prediction of the forces attri- 
butable to drag and to virtual mass is quite complex for flows in which 
separation takes place. At moments of maximum acceleration and zero 

velocity, there is no wake, and the pattern of flow closely resembles 

that for irrotational flow without separation. The foregoing dis- 

cussion and analysis explain the reason for the selection of a simple- 
harmonic wave motion in the determination of virtual mass of a circular 


cylinder or, for that matter, of any other object. 
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Ve Results. 

The following tables present representative values obtained during 
the course of the experiments. The folmentne procedure was used in the 
determination of the inertia coefficient: 
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Figure 5. Sample recording of amplitude of wave and 
force on cylinder. 


1. Inserting the measured amplitude of the wave and the appropri- 
ate values of k, H, and z into equation (33), noting that the 
amplitude was taken at its maximum when cos (oe t)= 1, the 
amplitude of the acceleration was calculated; 


2. Using the corresponding measured values of the amplitude of 
force, the inertia coefficient was calculated from equation 


(6). ; 
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TABLE I 
Results of system check with fresh water 
Cylinder diameter: 1 inch 
Water depth: 19 inches 
Height of cylinder above bottom: 4.375 inches 


RUN NUMBER ONE 


Point Force on Amplitude Inertia 
cylinder of wave Coefficient | 
wo a a a wo ounds inches 
i 0124 281 2.154 
2 00976 soo 1,906 
3 0092 2La7 2.054 
4, .008 » 2102 1.842 
5 OOH Z «lev 1.875 
6 .006 141 2.077 
7 006 »141 2.077 
8 0044 wL05 2.046 
9 0036 .0937 1,876 
10 00264 0625 2,062 
RUN NUMBER TWO 
1 .018 4225 2.080 
2 016 .406 1.924 
3 0145 . 360 1.966 
4 0128 . 328 1,905 
5 0116 202 2.008 
6 0104 a0 2.031 
7 0092 » 200 L.9L1 
8 .0083 .208 1.948 
9 00575 »141 1.991 
10 0048 .1158 2.023 
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TABLE I (continued) 
Results of system check with fresh water 


RUN NUMBER THREE 


Point Force on Amplitude Inertia 
cylinder of wave Coefficients 
(pounds ) (inches) 
1 0232 750 2.085 
2 029 wy 48, 1.968 
3 O28, . 666 1.93% 
4 0255 .610 2.043 
a 0280 ey) 1.93e 
6 « O2mi2 .506 2.047 
7 ,0206 . 500 2.015 
8 .0185 445 2.0g2 
9 Ode 406 2.046 
10 0 375 12955 
ll 0136 344 1.982 
IZ .0124 ole 1.942 
13 .0114 281 1.982 
14 0098 250 1.916 
15 0094 - on 2.098 
16 . 0084 203 2.022 
17 0078 « Lae 5 2.033 
18 0065 . L562 2.034 
19 0058 ~ 1407 2.015 
20 0052 - ize 2.033 
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TABLE II 
Results of system check with fresh water 
Cylinder diameter 1.5 inches 


Water depth: 19 inches 
Height of cylinder above bottom: 6.5 inches 


Point Force on Amplitude Inertia 

cylinder of wave Coeffictent 
(pounds) _ (inches) 

1 0146 ,150 1.957 

2 024 125 1.930 

3 026 0275 1,901 

4 057) . 388 1.943 

> .040 »413 1,947 

6 0415 »425 1,963 

7 020 » 200 2.010 

8 0218 feo 1.948 

9 0300 Be siya) 1,930 

10 034 .338 2.022 

TABLE ITIL 


Cylinder diameter: 2 inches 
Water depth: 18.875 inches 
Height of cylinder above bottom: 8.5 inches 


1 .058 en 2 2.097 
Z »055 272 1,988 
2 »052 »247 2.070 
b .030 »147 2.007 
5 029 » 136 2.097 
6 . 080 »358 2,198 
7 .070 » 321 2.144 
8 062 0297 2.053 
9 ,038 » 185 2.020 
10 024 2124 1.903 
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TABLE IV 
Results in stratified salt-water 


Density gradient: .000020 #/in’ 

Water depth: 19.5 inches 

Cylinder diameter: 1 inch 

Height of cylinder above bottom: 6.5 inches 


Point Force on Amplitude Inertia 
cylinder of wave Coefficient 
(pounds) (inches) 
li .0309 . Sil 2.83 
2 .0298 -48 2.90 
3 ~0270 -46 2.74 
4 .0258 425 2.84 
5 .024 -40 2.80 
6 .0222 Pe) 2.80 
7 .0205 335 2.74 
8 .0195 Re 236) 
9 » OMT 2 30 2.68 
10 .0120 . oe 2 167 
11 .0116 .20 2.70 
12 .010 « NS 2.67 
is 009 16 2.63 
14 008 «13 2.88 
1S .0073 . ae 2.84 
16 .032 «Ne 2.88 
17 .0289 495 2. tie 
18 026 43 2 oe 
19 20245 a5 2.89 
20 oONK?'S A) Z ae 
21 . 160 26 2. on 
22 .0142 ~ 24 2.00 
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TABLE V 
Results in stratified salt-water 


Density gradient: .000020 #/in’ 

Water depth: 19.5 inches 

Cylinder diameter: 1.5 inches 

Height of cylinder above bottom: 6.5 inches 


Point Force on Amplitude Inertia 
cylinder of wave Coefficient 
(pounds) (inches) 
1 0475 £38 2.64 
2 048 36 2.48 
3 046 36 2.68 
4 044 - 362 2.58 
5 035 26 2.60 
6 Oe22 28 2.68 
7 029 24 2 ae 
8 0263 20 2 eG 
o . 022 18 2.54 
10 .0198 . 16 2.508 
11 O75 461 2.56 
12 Pal 8 74 20 2.56 
13 068 254 2 4G2 
14 064 ow 2.66 
15 .0585 46 2.65 
16 055 44 2.60 
17 . O62 415 2. Gm 
18 049 40 2.56 
19 046 . 365 2.62 
20 -0865 29 2.62 
21 028 22 2.65 
22 . O23 18 2.66 
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TABLE VI 
Results in stratified salt-water 


Density gradient: .000020 #/in* 

Water depth: 19.5 inches 

Cylinder diameter: 2 inches 

Height of cylinder above bottom: 6.5 inches 


Point Force on Amplitude Inertia 
cylinder of wave Coefficient 
(pounds) (inches) 
1 096 44 2.56 
2 ,088 -40 2.39 
3 . 080 2 36 2.60 
4 073 ~345 2.48 
5 0645 29 2.60 
6 0565 265 2.90 
7 054 26 2.44 
8 046 ae 2.45 
9 0392 18 2.95 
10 036 « Le 2.64 
ll 0255 oo 2.46 
12 0295 135 2.56 
13 .0275 ~12 2.68 
14 025 ell 2.66 
15 .0225 .10 2.64 
16 » 146 66 2.90 
17 107 47 2.67 
18 .078 -J2 2.61 
i 0405 .18 2.63 
20 .0265 12 2.58 
21 094 42 2.62 
22 . 038 Be a) 2.62 
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TABLE VII1 
Results in stratified salt-water 


Density gradient: .000042 #/int 

Cylinder diameter: 1 inch 

Water depth: 19.5 inches 

Height of cylinder above bottom: 6.5 inches 


Point Force on Amplitude Inertia 
cylinder of wave Coefficient 





1 
2 .0133 21 2.85 
3 »0104 18 2,66 
4 O10 16 2,88 
5 .0093 Pa 2.85 
6 008 AF 2.83 
7 0435 ys! 2.82 
8 0405 67 2.78 
9 0355 .58 2.82 
10 0265 44 20/1 
11 026 42 2.85 
12 0245 40 2.82 
13 .0232 38 2,81 
14 022 36 2,81 
15 . 020 39 2,688 
16 919 30 2.91 
17 0172 28 2.83 
18 .0154 a 2.83 
19 O14 23 2,80 
20 .9132 622 2,76 





TABLE VIIL 
Results in stratified salt-water 


Density gradient: .000042 #/in’ 

Cylinder diameter: 1.5 inches 

Water depth: 19.5 inches 

Height of cylinder above bottom: 6.5 inches 


Point Force on Amplitude Inertia 
cylinder of wave Coefficient 
(pounds) (inches) 
1 06 047 2.60 
2 058 44 2.09 
3 055 ~43 2.61 
4 052 -40 2.65 
5 050 38 2.69 
6 046 » 36 2.61 
7 0435 34 2.61 
8 ~0417 ry 5: 2.58 
9 0403 305 2.69 
10 «0375 200 2o00 
11 0356 28 2.39 
12 .0325 myA®) 2.56 
13 .0295 24 2.51 
14 026 «aes 2.52 
15 .GZ> 20 Zao 
16 0158 P| 2.48 
LZ -0108 02 2.45 
18 084 -63 2 ne 
19 063 oO parse 
20 049 ae 2.56 
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TABLE IX 


Results in stratified salt«water 


Density gradient: .000042 #/in’ 

Cylinder diameter: 2 inches 

Water depth: 19.5 inches 

Height of cylinder above bottom: 6.5 inchés 








Point Force on Anplitiude Inertia 
cylinder of wave Coefficient 
jounds) _.._(ifte , em oe 
1 241 , 64 2.54 
2 130 .58 2,58 
3 ,110 +50 2:53 
4 . 103 345 2,63 
5 ,092 541 2.58 
6 .078 » 36 2,49 
7 . 069 6315 2,52 
8 , 060 27 2.56 
9 .052 245 2.44 
10 042 20 2.42 
11 ;0385 .18 2.46 
12 036 :17 2.44 
13 ;0315 .15 2.42 
14 100 ,45 2,56 
15 080 »365 2,52 
16 067 29 2.66 
17 .0565 26 2.50 
18 0485 22 2,535 
19 037 16 2,66 
20 026 ;12 2,49 
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TABLE X 


A comparison of theoretical and experimental inertia coefficients. 


Experimental values are means of numerous points taken during several 
runs. Note the very slight effect of the cylinder size and density 


stratification on the theoretical values of inertia coefficients. 


Inertia coefficients 


Cylinder size Theoretical Experimental 
(inches) value mean 
Fresh water 
1.0 2.000 1.998 
es 2.000 1.910 
2 0 2.000 2.094 


First stratification, gradient 0.00002 sesake 


Ie 2.500275 2576 
1.5 2.500413 2 «61 
2 10 2.500550 2.565 


Second stratification, gradient 0.000042 Bis il 


1.0 2.500715 2.834 
1.5 2.501074 2.508 
2.0 2.501433 2.509 
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VI. Discussion of Results and Conclusions. 

The experimental data tabulated in tables I-X have shown that the 
inertia coefficient for circular cylinders depends on the density of the 
fluid and to a lesser extent on the density gradient and the diameter of 
the cylinder. For reasons yet unknown the inertia coefficients for the 
one inch cylinder are about 15% higher than those for the 1.5 and 2 inch 
cylinders. The values calculated with the appropriate values for the 
density, density gradient, and the diameter of the cylinder through the 
use of equation (24) gave an inertia coefficient of 2.5 as shown in 
Table X. This value varied very slightly with either the density gradi- 
ent or the cylinder size. It is, therefore, apparent from a comparison 
of the theoretical value with those obtained experimentally that the 
theoretical value is somewhat lower than those predicted experimentally. 
Among several reasons which could be provided, one is that the stratified 
fluid medium is confined and that this confinement, unlike in the case 
of homogeneous fluids, affects more adversely the blockage effect cited 
earlier. Additional errors may have come from the deviations of the 
density gradient from uniformity and from the presence of the connecting 
arms extending through a stratified media. Nevertheless, the agree- 
ment between the theory and experiments is sufficiently good to conclude 
that the inertia coefficient for a circular cylinder immersed in a strati- 
fied fluid media is about 25% larger than that in a homogeneous media. 
The consequences of the increase of the inertia coefficient are rather 
apparent when one considers the fact that in such media as the oceans, 
and in particular in layers close to the thermocline, the acceleration 


of a body requires a proportionally larger amount of power. 
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In the foregoing we have discussed so far the theoretical and 
experimental values of the inertia coefficients as they were calcu- 
lated or measured, but did not elaborate on the reasons as to why such 
an increase should occur. A careful analysis of the flow about a circu- 
lar cylinder immersed in a stratified fluid medium (See Fig. 6) shows 


that there are six instead of two stagnation points about the cylinder. 





Figure 6. Flow of a stratified fluid about a circular cylinder. 
Furthermore, there are on both sides of the cylinder two roughly tri- 
angular eddies symmetrically located with respect to either axis. It 
is apparent that the presence of additional stagnation points and the 
wedge shaped fluid body in the front and rear of the cylinder would 
necessarily increase the added mass and hence the inertia coefficient. 

In conclusion, it suffices to state that the inertia coefficient 
for cylinders immersed in a stratified fluid medium depends on the 
size of the cylinder, the density gradient, and the density of the 
fluid, and that the inertia coefficient is about 25% larger than that. 


of a cylinder immersed in an homogeneous medium. 
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VII. Recommendations for further Studies. 

It ig strongly recommended that the project deaeribed herein be 
extended to flat plates, spheres, and other three dimeneional bodiea 
through the use of the same experimental apparatus. Even though an analy= 
sis of the stratified flow about a three dimenaional bady may be rather 
difficult, it is believed that an analysis similar to that presented 
here may be carried out for a flat plate vertically placed in the reger- 
voir, 1.¢., normal to the direction of the horizontal component ef the 


velocity. 
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